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ABSTRACT 

We determine the distances to the z ~ 0.55 galaxy clusters MS 0451.6 — 0305 and CI 0016 + 16 from a 
maximum likelihood joint fit to interferometric Sunyaev-Zel'dovich effect (SZE) and X-ray observations. 
We model the intracluster medium (ICM) using a spherical isothermal j3 model. We quantify the 
statistical and systematic uncertainties inherent to these direct distance measurements, and we determine 
constraints on the Hubble parameter for three different cosmologies. For an Um — 0.3, ft a = 0.7 



cosmology, these distances imply a Hubble constant of 63 



+ 12 +21 

-21 



km s 1 Mpc 1 , where the uncertainties 



correspond to statistical followed by systematic at 68% confidence. The best fit H is 57 km s 1 Mpc 1 
for an open Qm — 0.3 universe and 52 km s _1 Mpc -1 for a flat fijyj — 1 universe. 

Subject headings: cosmic microwave background — cosmology: observations — distance scale — 
galaxies: clusters: individual (MS 0451.6- 0305; CI 0016+16) — techniques: 
interferometric 



1. INTRODUCTION 

Analysis of Sunyaev-Zel'dovich effect (SZE) and X-ray 
data from a cluster of galaxies provides information that 
can be used to determine the diatancc to the cluster, in — 



recent success of SZE observations is due to advances in 
instrumentation and observational strategy. Recent high 
signal-to-noise ratio detections have been made with sin- 



dopondont of tho oxtragalaotio diutc 



In tho 



gle dish observations at radio wavelengths (Birkinshaw & 
Hughes 1994]; iHerbiget al. 1995| ; pVlyers et al. 1997j;|Hu°lic " 



oar h i 



rontioo, Si 



and Zol'doTfioh (1070, 1070) Dug 

goetod that cosmic microwavo background (CMB) photonc 



fc Birkinshaw 1998 ), millimeter wavelengt hs"!' Holzapfe 



inverse Compton scattering off the electrons in the hot 
{f < 10 koV) intracluster medium (ICM) trapped in tho 



et al. 1997a. b|; Pointccoutcau ct al. 1999[) and submil- 
limct cr wavelengths ( Lamarrc ct al. 1998| ; [komatsu et al 



1999). Interferometric ob servations have produced high 



potential well of the cluster, would cause a small (<; 1 



quali t y images of the SZE ( Jones et al. 1993 [_ Grainge et al 
1993|; barlstrom et al. 1996j; Parlstrom et al. 1998j ; |Saun- 



111K) dLMuiLiuii ill LhL OMD ipLUmiii, mm kiiuwn as Qil ders et al. 19991; brainge et al 



aunyacv-Zcl'dovich effect (SZE). The distortion appears 

as a decrement for frequencies <; 218 GHz (A J> 1.4 mm) 



1999). On the theoretica 



land 



218 GHz. The SZE 



as ail increment, tor frequencies ;> 
' signal is proportional to the pressure integrated along the 

line of sight through the cluster, AT ~ J n e T e d£, where n e 
is the electron density of the ICM and T e is the electron 
temperature. The X-ray surface brightness can be written 
as S x ~ J n 2 e K eli d£ where A eH is the X-ray cooling func- 
tion, which depends on temperature and metallicity. It 
was soon realized that one can determine the distance to 
the cluster by capitalizing on the different dependencies 
on density, n e , with some assumptions about the geom- 
etry of the cluster. This is a direct distance based only 
on relatively simple cluster physics and not requiring any 
standard candles or rulers. 

The SZE signal is weak and difficult to detect. The 



side, there is substanti al literature on relativistic correc- 
tions to both the SZE (|Rcphacli fe Yankovitch 1997j lltoh 



et al. 1998; Challinor fc Lasenby 199S ; Nozawa et al. 1998a 



Sazonov fc Sunyaev 1998|) and X-ray brcmsstrahlung 



(IHughes fc Birkinshaw 1998| : |Rcphaeli fc Yankovitch 1997 



Nozawa et al. 1998 q] ). Todatc, there are about a dozen 



estimates of H Q based on combining X-ray and SZE data 
for individual clusters (see Birkinshaw 1999 for a compre- 
hensive review). 

We present a new analysis, wherein we perform a joint 
maximum-likelihood fit to both interferometric SZE and 
X-ray data. This method takes advantage of the unique 
properties of interferometric SZE data, utilizing all the 
available image data on the ICM. This is the first time 
SZE and X-ray data have been analyzed jointly. We ap- 
ply this method to observations of MS 0451.6 — 0305 and 
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Table 1 
Radio Point Sources 





RA 


DEC 


^28.5 


^15 


F 5 




Field 


(J2000) 


(J2000) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


MS 0451 


04 h 54 m 22 s 


-03° 01' 26" 


1.88 






14.9 a 


CI 0016 


00 h 18 m 31 s 


+16° 20' 45" 


9.07 


25.0 b 


84.5 b 


267 b 



a From Condon et al. (1998). 
b From Moffet & Birkinshaw (1989). 



CI 00 16 + 16, massive clusters at rcdshift z = 0.55 (Don 
ahue k Stockc 19951; ICarlbcrg et al. 1994|) and z = 5455 



( Neumann fc Bohringcr 1997| ; Dressier fc Gunn 1992 ), re- 
spectively^ 



in 1997 between June 21 and July 22 for 8 hours with nine 
6.1 m telescopes, both years with an 800 MHz bandwidth 
centered at 28.5 GHz. 



We describe the data and reduction in jj |2|, the analysis 



Th e data are reduced using the MIRIAD ( Bault et al 
1995) software package at BIMA and using MMA (| cov- 



atic uncertainties in § |]. Section contains a discussion 
of the results and future prospects. All uncertainties are 
68.3% confidence unless explicitly stated otherwise. 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. Interferometric SZE Observations 

The extremely low systematics of interferometers and 
their two-dimensional imaging capability make them well 
suited to study the weak (<; 1 mK) SZE signal in galaxy 
clusters. Over the past several summers, we outfitted the 
Berkeley Illinois Maryland Association (BIMA) millime- 
ter array in Hat Creek, California, and the Owens Valley 
Radio Observatory (OVRO) millimeter array in Big Pine, 
California, with centimeter wavelength receivers. Our re- 
ceivers use cooled (~ 10 K ) High Electron Mobility Tran- 
sistor (HEMT) amplifiers ( Pospieszalski et al. 1995| ) oper- 
ating over 26-36 GHz with characteristic receiver temper- 
atures of T rx ^11-20 K, over 28-30 GHz, the band used 
for the observations presented here. When combined with 
the BIMA or OVRO systems, these receivers obtain typ- 
ical system temperatures scaled to above the atmosphere 
of T sys ~ 45 K and as low as 34 K. Most telescopes are 
placed close together in a compact configuration to probe 
the angular scales subtended by distant clusters (~ 1'), 
but telescopes are always placed at longer baselines for 
simultaneous detection of point sources. Every half hour 
we observe a radio point source, commonly called a phase 
calibrator, to monitor the system gains for about two min- 
utes. 

MS 0451 was observed at OVRO in 1996 during May 
and June for 30 hours with six 10.4 m telescopes using 
two 1 GHz channels centered at 28.5 GHz and 30.0 GHz 
(2 GHz bandwidth). CI 0016 was observed at OVRO in 
1994 between June 16 and July 4 for 87 hours with five 
10.4 m telescopes and a 1 GHz bandwidth centered at 
28.7 GHz and in 1995 between July 24 and July 28 for 13 
hours using five 10.4 m telescopes and two 1 GHz chan- 
nels centered at 28.5 GHz and 30.0 GHz. CI 0016 was 
also observed at BIMA in 1996 between September 6 and 
September 18 for 29 hours with six 6.1 m telescopes and 



methc d"*in §|3|, and present the results and possible system- illc et al. 199*3 ) at OVRO. In both cases, data are removed 



when one telescope shadows another, when cluster data 
are not straddled by two phase calibrators, when there 
are anomalous changes in instrumental response between 
calibrator observations, or when there is spurious correla- 
tion. For absolute flux calibration, we use observations of 
Mars, on the assumption that we know its true brightness 
temperature from the Rudy (1987) Mars model. For ob- 
servations not containing Mars, calibrators in those fields 
are bootstra pped back to the nearest Mars calibration (see 
Grego 1999 for more details). The observations of the 



phase calibrators over each summer give us a summer- 
long calibration of the gains of the BIMA and OVRO in- 
terferometers. They both show very little gain variation, 
changing by less than 1% over a many-hour track, and the 
average gains remain stable from day to day. 

An interferometer samples the Fourier transform of the 
sky brightness rather than the direct image of the sky. The 
final products from the interferometer are the amplitudes 
of the real and imaginary components of the Fourier trans- 
form of the cluster SZE distribution on the sky multiplied 
by the primary beam of the telescope. The SZE data files 
include the positions in the Fourier domain, which depend 
on the arrangement of the telescopes in the array, the real 
and imaginary components, and a measure of the noise in 
the real and imaginary components. The Fourier conju- 
gate variables to right ascension and declination are com- 
monly called u and v, respectively, and the Fourier domain 
is commonly referred to as the u-v plane. 

The finite size of each telescope dish imposes an almost 
Gaussian attenuation across the field of view, known as the 
primary beam. The primary beams are constructed from 
holography data taken at each array. The main lobe of the 
primary beams can be approximated as Gaussian with a 
full width at half maximum (FWHM) of 4'. 2 for OVRO 
and 6'. 6 for BIMA at 28.5 GHz. We use the primary beam 
profiles made from the holography data for our analysis. 

The primary beam sets the field of view. The effec- 
tive resolution, called the synthesized beam, depends on 
the sampling of the u-v plane of the observation and is 
therefore a function of the configuration of the telescopes. 



3 




4 h 54 m 20 s 16 s 12 s 8 s 4 s h 18 m 48 s 42 s 36 s 30 s 24 s 

Right Ascension (J2000) Right Ascension (J2000) 

Fig. 1. — SZE (contours) and X-ray (grey scale) image overlay for MS 0451.6 — 0305 and CI 0016 + 16. The SZE images have an rms of ~ 30 
fiK (MS 0451) and ~ 46 fiK (CI 0016). The contours are multiples of 2 a and negative contours are shown as solid lines. The FWHM ellipse 
of the synthesized beam is shown in the lower left corner of each panel. The different sizes of the SZE images are due to the smaller dishes 
of the BIMA telescopes (used for CI 0016) which can be more closely packed and thus sample smaller u-v radii than the OVRO telescopes 
(used for MS 0451). The X-ray images are PSPC raw counts smoothed with a Gaussian with a = 15". The peak for the MS 0451 image is 
23 counts and 50 counts for the CI 0016 image. 



The c [uster SZE signal is largest on the shortest baselines 
(largest angular scales). The shortest possible baseline is 
set by the diameter of the telescopes, D. Thus we are not 
sensitive to angular scales larger than about X/2D, which 
is ~ 2'. 8 for BIMA observations and ~ T.7 for OVRO 
observations. The compact configuration used for our ob- 
servations yields significant SZE signal at these angular 
scales, but the interferometer is not sensitive to larger an- 
gular scales. Because of this spatial filtering by the inter- 
ferometer, it is necessary to fit models directly to the data 
in the u-v plane, rather than to the deconvolved image. 

Point sources are identified from SZE images created 
with DIFMAP ( Pearson et al. 1994 ) using only the long- 
baseline data (;> 2000 A) and natural weighting. Approx- 
imate positions and fluxes for each point source are ob- 
tained from this image and used as inputs for the model 



et al. 1998). The point source in the CI 0016 field is 339" 



fitting discussed in § |3.2| . The data are separated by obser- 



vatory, frequency, and by year to allow for temporal and 
spectral variability of the point source flux. One point 
source is found in both the MS 0451 and CI 0016 fields. 
The point source pos itio ns and fluxes are from the model 
fitting described in § 3^ and summarized in Table [l]. Our 
positions agree very well with the NRAO VLA Sky Survey 
(NVS S) source positions (Condon et al. 1998 ) 



from the pointing center and is only seen in the BIMA 
data since the OVRO primary beam attenuation places it 
out of the OVRO field of view. The flux of this source is 



measured to be I.OI+oJq mJ y at 28 - 5 GHz from the 1997 
BIMA data, which when corrected for the primary beam 
attenuation, is an intrinsic flux of 9.07^'gQ. This source 
corresponds to source 15 from the survey of this field done 
by Moffet & Birkinshaw (1989). They found this source to 
be 267 ± 3 mJ y at 1.44 GHz (264 mJy in the more recent 



NVSS survey; |Condon et al. 1998[ ), 84.5 ±1.1 mJy at 4.86 



GHz, and 25.0 ± 1.5 mJy at 14.94 GHz. We do not see the 
other two sources of Moffet & Birkinshaw (1989) within 
360" of the pointing center, sources 10 and 14 which have 

15 GHz fluxes of 0.56 and < 2.7 mJy, respectively. We do 
not present a flux for the point source in the 1996 BIMA CI 

0016 data because of a problem with the absolute calibra- 
tion of the array during observations early that summer. 
Though the overall normalization is uncertain, the data 
still provide shape information about the cluster. Table ^ 
summarizes the positions of the point sources and their 
fluxes at various frequencies. 

Figure [l] shows the SZE image contours overlai d on the 
X-ray imag es of these clusters. We use DIFMAP (Pearson 



Th( | MS 0451 field point source is located 172" from the et al. 1994| ) to produce the naturally weighted SZE im 



pointing center with a measured flux of 0.50^^ mJy at 
28.5 GHz. Correcting for the primary beam attenuation 
appropriate for this offset, the intrinsic point source flux 
is 1.88lQ2g m Jy In the 30 GHz channel, this point source 
has a flux of 0.36^0 07 m Jyj an d l-35ip'2g after correcting 
for the primary beam. This point source was foun d in the 
NVSS survey with a flux of 14.9 mJy at 1.4 GHz ( pondorj 



ages. The point sources are subtracted from the data and 
a Gaussian taper applied to emphasize brightness vari- 
ations on cluster scales before the image is deconvolved 
(CLEANed). For the MS 0451 OVRO data, we apply 
a 1200 A half-power radius Gaussian taper before decon- 
volving the image. This yields an elliptical Gaussian fit 
of 48" x 70" for the synthesized beam (effective resolu- 



4 



tion) and a rms of ~ 68 fiJy beam , corresponding to a 
Rayleigh- Jeans (RJ) brightness sensitivity of ~ 30 fiK. For 
CI 0016 we use the 1996 and 1997 BIMA data with a 1000 
A half-power radius Gaussian taper giving a 81" x 101" 
synthesized beam and a rms of ~ 250 /xJy beam" 1 , cor- 
responding to a ~ 46 fiK RJ brightness sensitivity. The 
SZE image contours are multiples of twice the rms level for 
each image. Images made with the 1994 OVRO CI 0016 
data were presented in Carlstrom, Joy, & Grego (1996). 

We stress that these images are made to demonstrate 
the data quality. The actual analysis is done in the Fourier 
plane, where the noise characteristics of the data and the 
spatial filtering of the interferometer arc well understood. 
The SZE and X-ray image overlays in Figure |l| show that 
the region of the cluster sampled by the interferometric 
SZE observations and the X-ray observations is similar. 

2.2. X-ray Observations 

We use archival Rontgen Satellite (ROSAT) data from 
both the Position Sensitive Proportional Counter (PSPC) 
and High-Resolution Imager (HRI) instruments. MS 0451 
was observed with the PSPC in 1993 over March 5-7 for 
15,439 s of live-time and by the HRI in 1995 over Septem- 
ber 3-19 for 45,864 s of live-time. There are approximately 
1200 cluster photons collected in both the PSPC and HRI 
observations. CI 0016 was observed with the PSPC in 
1992 over July 11-18 for a live-time of 41,589 s and by the 
HRI in 1995 between June 17 and July 5 for a live-time 
of 70,228 s. The PSPC data contains about 3200 cluster 
photons and the HRI has about 1500. 

We use the Snowden Extended Source Analysi s Soft- 
ware (ESAS) ( |5nowdcn ct al. 1994| ; |Snowdcn 199g| ) to re- 
duce the data. We use this software to generate a raw 
counts image, a noncosmic background image, and an ex- 
posure map for the HRI (0.1-2.4 keV) data and for each 
of the Snowden bands R4-R7 (PI channels 52 — 201; ap- 
proximately 0.5 — 2.0 keV) for the PSPC data, using a 
master veto rate (a measure of the cosmic-ray and 7-ray 
backgrounds) of 200 counts s _1 for the PSPC data. We 
examine the light curve data of both instruments looking 
for periodic, anomalously high count rates (short-term en- 
hancements) and for periods of high scattered solar X-ray 
contamination. None are found. The Snowden software 
produces 512 x 512 pixel images with 14.947" pixels for 
the PSPC and 5.0" pixels for the HRI. For the PSPC, 
final images for all of the R4-R7 bands together are gener- 
ated by adding the raw counts images and the background 
images. Each Snowden band has a slightly different effec- 
tive exposure map and there is an energy dependence in 
the point spread function (PSF). Thus, we generate a sin- 
gle exposure image and a single PSF image by combining 
cluster photon-weighted averages of the four exposure im- 
ages and the four PROS ( Worrall et al. 1992; Conroy et al 



1993) generated on-axis PSF images. The cluster photon- 
iting is determined using the background subtracted 
detected photons within a circular region centered on the 
cluster. The region selected to construct the weights is 
the largest circular region encompassing the cluster which 
contains no bright point sources. For MS 0451, we use a 
12 pixel radius and for CI 0016 we use a 15 pixel radius. 
X-ray images with SZE image overlays of MS 0451 and 



CI 0016 are shown in Figure [l[ The gray-scale images are 
the PSPC "raw" counts images smoothed with a Gaussian 
with ct = 15". The peaks of the images are 23 counts (MS 
0451) and 50 counts (CI 0016). 

For MS 0451, we use the emission- weighted temper- 
ature, galactic absorption, and metallicity from Don- 
ahue (1996). She found a best-fit X-ray temperature of 
T e = 10.4l g keV with a galactic absorption column 
density of N H — 3.0^o'3 x 10 20 cm" 2 and a metallicity 

of 0.15^0 07 s °l ar implied by the iron abundance from a 
joint analysis of ASCA and PSPC data. This temper- 
ature is consistent with the Mushotzky & Scharf (1997) 

value of T e = 10.17t°:?6 keV - For C1 0016 > we a^P* 
the Hughes and Birkinshaw (1998) results. They found 
T P = 7.55^0 58 keV with a galactic absorption column den- 



sity of N H = 5.591^36 x 10 2U cm" 2 and a metallicity of 
0-07to.l\ solar from a joint analysis of ASCA and PSPC 
data. This temperature agrees with a more recent analy- 
sis by Furuzawa et al. (1998) who found T e = 8-0±g | keV. 
Unlike the Hughes & Birkinshaw analysis, this analysis 
did not include the PSPC data which is sensitive to the 
column density. 

2.2.1. X-ray Cooling Function 

The Raymond-Smith (1977) code calculates the 
electron-ion bremsstrahlung contribution in the non- 
relativistic limit using the Gaunt factors of Karz as fc Lat- 
ter (1961). Recent l y it has been pointed out ( Rephacli 



fc Yankovitch 1997| ftughes fc Birkinshaw 1998j ) that rel- 
ativistic corrections are important for a precise deter- 
mination of angular diameter distances, and therefore 
also of H Q . Though dubbed "relativistic" corrections, 
the corrections go beyond only relativistic effects and in- 
clude 1) relativistic corrections to the electron distribu- 
tion function, 2) relativistic and spin corrections to the 
non-relativistic electron-ion bremsstrahlung cross section, 
3) inclusion of electron-electron bremsstrahlung (all three 
of relative order k B T e /m e c 2 ), and 4) first-order Born ap- 
proximation corrections to electron-ion bremsstrahlung (of 
order (Ry /k B T e y^ 2 , where Ry = 13.6 eV is the ioniza- 
tion energy of hydrogen). When applied, these corrections 
provide a better than 1% acc urate calculation of thermal 
bremsstrahlung (Gould 1980). Gould provides results for 
both the integrated energy-loss rate and the spectral cool- 
ing function from thermal bremsstrahlung. 

Following Hughes & Birkinshaw (1998), we use the cor- 
rections to the spectral cooling function rather than the to- 
tal energy-loss rate used by Rephaeli & Yankovitch (1997) 
because the calculated cooling function comes from inte- 
grating the spectral cooling function over the fairly nar- 
row ROSAT energy band redshifted to the cluster frame. 
These corrections may also affect the T e derived from fits 
to X-ray spectral data because these corrections modify 
the shape of the X-ray spectrum. Hughes & Birkinshaw 
(1998) found a ~ 1% change in the best fit T e for the Coma 
cluster (T e ~ 8 keV) when applying these corrections to 
their spectral fits. As a check we also verified that the spec- 
tral cooling function formula from Gould when integrated 
over energy agrees with the total energy- loss result 10 . 
To calculate the X-ray spectral cooling function, we use 



°We have verified the misprint in Gould (1980) discussed by Hughes & Birkinshaw (1998). This misprint combined with the slightly different 
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Table 2 

X-ray Cooling Functions for the PSPC 



Cluster 


AeHO 

(erg s _1 cm 3 ) 


A dot 

J1 eHo 

(cnt s" 1 cm 5 ) 


£ 

,erg s _1 cm -2 
1 cnt s" 1 > 


MS 0451 
CI 0016 


6.95 x 10~ 24 
6.91 x 10~ 24 


3.26 x 10~ 13 
3.00 x 10~ 13 


1.37 x 10" 11 
1.49 x 10~ n 



a Rajmond-Smith (1993 Sep 21 version) thermal plasma 
mode l with its brcmsstrahhmg component replaced with 



et al. 1999a) 



Goul d/s bremsstrahlung calculation including the correc- 
tions discussed above. We use the Anders & Grevesse 
(1989) meteoritic abundances as the solar values, scaling 
the abundances of elements heavier than He by the metal- 
licity of the cluster. We calculate the absorption from cold 
Galactic gas using the photoelectric cross sections from 
Balucinska- Church & McCammon (1992) (including the 
updated He absorption; 1993 Sep 23 version) for Anders 
& Grevesse solar abundances. 

We integrate the modified Raymond-Smith spectral 
model over the redshifted ROSAThand (0.5-2.0 keV in the 
detector frame) to determine the cooling function in cgs 
units, A eHO - We also calculate the cooling function in de- 
tector units by multiplying the modified Raymond-Smith 
spectrum by the response 11 (includes effective area and 
energy resolution) of the instrument, dividing by the en- 
ergy of the photons (to convert to counts), and integrating 
to find the total cooling function, A^ . Comparing these 
two yields the detector to cgs unit conversion, £, after 
correcting by (1 + z) due to the difference between instru- 
ment counts and energy (S = A eH o/Ag^ /(l + z)). The 
cooling function results for the PSPC are summarized in 
Table @. The cooling functions are 1.052 and 1.046 times 
the Raymond-Smith "uncorrected" value for MS 0451 and 
CI 0016, respectively. 

3. METHOD 



The spherical isother mal (3 model is given by (Cavaliere 



fc Fusco-Fcmiano 1976| , 1978) 



3.1. Angular Diameter Distance Calculation 



The 
cluste 
mode 



calculation begins by constructing a model for the 
gas disLiibuLiuii. We use a spheiical Lsulheiiiial ft 
to describe the ICM. Within this context the clua - 



ter's characteristic scale along the line of sight is the same 
as the scale in the plane of the sky. This model is clearly in- 
valid in the presence of cluster asphcricities. Thus cluster 
geometry introduces an important uncertainty in SZE and 
X-ray derived distances. In general, clusters are dynam- 
ically young, are aspherical, and rarely exhibi t projected 
gas d istributions which are circular on the sky ( Mohr et al] 
1995). We currently can not disentangle the complicated 
cluster structure and projection effects, but numerical sim- 
ulations provide a good base for understanding these dif- 
ficulties. The effects of asphericity contribute significantly 
to the distance uncertainty for each cluster, but do not 
result in any significant bias in the Hubble parameter de- 



rived from a large sample of clusters ( [gulkancn 1999 : Mohi 



S (r) 



1 + — T 



-3/3/2 



(1) 



where n e is the electron number density, r is the radius 
from the center of the cluster, r c is the core radius of the 
ICM, and j3 is the power law index. With this model, the 
SZE signal is 



AT = f( x )T CMB D A / d( a T n e 



AT 1 + — 



(l-3/3)/2 



(2) 



where 



IS 



AT 
-4)(1 



the 

I- <$sz 



SZE decrement /increment, = 
5) (f(x) ~ * — 2 in the non-relativistic 



and Rayleigh- Jeans limits) is the frequency dependence of 
the SZE with x — hv /fcT CMB , S sz - E (x,T e ) is the relativistic 
correction to the f requency dependence, T CMB (=2.728 K; 
Fixsen et al. 1996| ) is the temperature of the CMB radia- 
tion, fc B is the Holtzmann constant, o~ T is the Thompson 
cross section, m e is the mass of the electron, c is the speed 
of light, AT is the central SZE decrement /increment, 9 
is the angular radius in the plane of the sky and 9 C the 
corresponding angular core radius, and the integration is 
along the line of sight £ = D A (. We apply th e relativis 



tic co rrections 5 SZE to fifth order in kT e /m e c ( [toh et al 



1998). The Itoh et al. results agree w ith other work 



bins 1997; Challinor fc Lascnby 1998| ) to third order where 
they stop. This correction decreases the magnitude of /( x ) 
by 3.7% for MS 0451 and 2.7% for CI 0016. The correction 
is slightly higher for MS 0451, as expected, because of its 
higher electron temperature. 
The X-ray surface brightness is 



1 



4it{1 + z) 4 



D A al(n e n H A 

e 



2X (l~6/3)/2 



e?. 



(3) 



where 
(erg s 
ter, n 



c is 
cm 

is the 



the 
2 



X-ray surface brightness in cgs units 
2 ), z is the redshift of the clus- 
density of the ICM, 



elemental compositions of the gas considered in Rybicki & 
between their total corrected energy-l oss rates (RL pg. 165 

11 Response matrices obtained from r tp://lcgacy.gsfc. nasa.gov/caldb/data/ 



arcmin J, z 
hydrogen number 

Lightman (1979; RL) and Rephaeli & Yankovitch (1997) explains the difference 
The RL value is correct for the pure hydrogen gas they consider. 

part of the HEASARC calibration database. 
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A eH = A eH (T e , abundance) is the X-ray cooling function 
of the ICM in the cluster rest frame in cgs units (erg cm 3 
s _1 ) integrated over the redshifted ROSATb&nd, and S xo 
is the X-ray surface brightness in cgs units at the cen- 
ter of the cluster. Since the X-ray observations are in 
instrument counts, we also need the conversion facto r be- 
tween detector counts and cgs units, E, discussed in § 2.2.1 
(S X o = S^H). The normalizations, AT and S xo , used in 
the fit include all of the physical parameters and geometric 
terms that come from the integration of the f3 model along 
the line of sight. 

One can solve for the angular diameter distance by elimi- 
nating n eo (noting that n K — n e a e / /i H where rij = p/ (ijm p 
for species j) yielding 



A = 



(ATp) 2 / m e C 2 \ A elio U e /u n 

S xo \k B T eo ) 4^/2/2 T 2 mb(7 2 (1 



(x) 



r(|/3-|) 



r(3/3 



r(3/3) 



(4) 



where T(x) is the Gamma function. Similarly, one can 
eliminate D A instead and solve for the central density n eo . 

3.2. Joint SZE and X-ray Model Fitting 

The SZE and X-ray emission both depend on the prop- 
erties of the ICM, so a joint fit to all the available data 
provides the best constraints on those properties. We per- 
form a joint fit to the interferometric SZE data and the 
PSPC and HRI X-ray data. Each data set is assigned a 
collection of parameterized models. Typically, SZE data 
sets are assigned a /3 model and point sources and X-ray 
images are assigned a (3 model and a cosmic X-ray back- 
ground model. This set of models is combined for each 
data set to create a composite model which is then com- 
pared to the data. 

Model parameters can be fixed, free to find their opti- 
mized values, or gridded. They can also be linked, forced 
to vary together among the data sets. In practice, 9 C and f3 
are linked between all data sets (both SZE and X-ray) and 
the central decrements AT are linked between the SZE 
data sets which are separated by season and array. We use 
a downhill simplex to sea rch parameter spa ce and maxi- 
mize the joint likelihood ( Press et al. 1992] ). The cluster 



position, /3, 9 C , S xo , AT , a constant cosmic background, 
radio point source positions, and point source fluxes are 
all allowed to vary. 

Each data set is independent, and likelihoods from each 
data set can simply be multiplied together to construct the 
joint likelihood. Likelihood ratio tests can then be per- 
formed to get confidence regions or compare two models. 
Rather than working directly with likelihoods, C, we work 
with S = — 21n(£). We then construct a A% 2 -like statistic 
from the log likelihoods, AS = S n — S re f where S re f is the 
reference S statistic, typically chosen to be the minimum 
of the S function, and S n is the S statistic where n pa- 
rameters differ from the parameters at the reference. The 
statistic AS is sometimes referred to as the Cash (1979) 
statistic and tends to a y 2 distrib ution with n degrees of 
freedom ( [Kendall fc Stuart 1979| for example). This AS 
statistic is equivalent to the likelihood ratio test and is 
used to generate confidence regions and confidence inter- 



vals with S re f = S m i n . For one interesting parameter, the 
68.3% (~ la) confidence level corresponds to A5" = 1.0. 

Because we arc interested only in differences in log likeli- 
hoods, AS, the model independent terms in the likelihoods 
are dropped. The log likelihoods are then 



"5 (A*? 



A/ 2 ) Wi for SZE data (Gaussian), (5) 



]Di\n(Mi) - Mi 



for X-ray data (Poisson), (6) 



where ARi and Ali are the differences between the model 
and data at each point i in the Fourier plane for the real 
and imaginary components respectively, Wi = l/of is a 
measure of the noise (Gaussian) of the real and imagi- 



nary components discussed in §2.1, and Mi and Di are 



the model prediction and data in pixel i. 

The interferometric SZE observations provide con- 
straints in the Fourier (u-v) plane, so we perform our 
model fitting in the u-v plane, where the noise properties of 
the data and the spatial filtering of the interferometer are 
well defined. The SZE composite model for both MS 0451 
and CI 0016 consists of a j3 model and a point source. The 
f3 model is computed in a regular grid in image space, mul- 
tiplied by the primary beam determined from holography 
measurements, and fast Fourier transformed to produce 
the u-v plane model. It is then interpolated to the u-v 
position for each data point. Point sources are computed 
analytically at each u-v data point of the observation and 
added to the (3 model in the u-v plane to construct the 
composite SZE model. The Gaussian likelihood (eq. ||) 
is calculated using the composite SZE model and the SZE 
data. During the fitting, the cluster center, 9 C , /3, AT , 
the point source positions, and the point source fluxes are 
all allowed to vary. 

Because the SZE is frequency dependent, a minor addi- 
tional detail comes from observations of a cluster at multi- 
ple frequencies, 28.5 GHz and 30 GHz. We input a central 
decrement appropriate for 30 GHz into the fitting routine, 
which then corrects the model for the actual frequency 
of the observation. The likelihood is calculated with the 
model appropriate for the observing frequency. This al- 
lows us to link the central decrement across data sets with 
different observing frequencies. 

The model for each X-ray data image includes a spheri- 
cal isothermal /? model plus a constant cosmic background. 
The model is convolved with the appropriate PSF, mul- 
tiplied with the exposure map, and then the noncosmic 
background is added pixel by pixel. Point sources are 
masked out. The logarithm of the Poisson likelihood 
(eq. ||) is then calculated. During the fitting, the clus- 
ter center, 9 C , (3, S xo , and the cosmic background are all 
allowed to vary. The PSF is generated by PROS and the 
exposure map and noncosmic background maps are those 
generated by the Snowden ESAS software discussed in 
§ |2^. The PSF and exposure maps for the PSPC Snowden 
bands R4-R7 a re c ombined in a cluster photon-weighted 
average (see § 2.2). Point sources are found using the 
ESAS detection algorithm with a 3 a detection criterion 
and masked out. Circular regions of typically 3 pixel ra- 
dius are placed on each point source and excluded from the 
calculation of the likelihood. These regions correspond to 
radii of - 45" for the PSPC and 15" for the HRI. As 




Radius (arcmin) 



Radius (arcmin) 



Fig. 2. — Comparison of the PSPC radially averaged surface brightness profile (points with error bars) with the best fitting isothermal 
spherical /3-model plus background (histogram) for a) MS 0451 and b) CI 0016. The lower panel shows the residuals in units of the standard 
deviation. The best fit model is a good fit to the data over the entire range of radii considered in the fit. There is no evidence of excess 
emission near the core for either cluster, the signature of cooling flows. 



Table 3 
ICM Parameters 



Cluster 




0c 

(arcsec) 


Cdct 

xo 

(cnt s^ 1 arcmin- 2 ) 


(erg s _1 cm -2 arcmin -2 ) 


AT 
( M K) 


MS 0451 
CI 0016 


0.806±S5g 
0.7491^ 


34.7 + f-§ 
42 3+ 2 4 


6.96j£S? xlO- 2 

4.i4ig;" xio- 2 


9.56+git xlO- 13 

Q.ntiH xio- 13 


-14311 H 
-1242±ig| 



a check, the image of the cluster excluding the masked 
regions is visually inspected. Increasing the size of the 
masked regions does not significantly alter the best fit pa- 
rameters, including the cosmic background. For the model 
fitting we use a region centered on the cluster with a 64 
pixel radius, corresponding to ~ 16' for the PSPC and 
~ 5' for the HRI. Using a larger fitting region does not 
change the best fit model parameters significantly. 

When allowed to vary separately, the best-fit central 
surface brightnesses for the PSPC and HRI are consistent 
within their uncertainties when compared in cgs units. We 
linked the central surface brightnesses between the PSPC 
and HRI in cgs units, using E to convert to counts before 
comparing with the X-ray images. The linked S X o case 
gives an insignificant change in the S statistic compared 
to the case where the PSPC and HRI normalizations are 
allowed to vary individually and removes one free param- 
eter. The consistent central surface brightnesses from the 
PSPC and HRI observations present an interesting test of 
the relative calibration of the two instruments. 

4. DIRECT DISTANCES AND THE HUBBLE CONSTANT 

The results from our maximum likelihood joint fit to the 
SZE and X-ray data are summarized in Table || for both 
clusters. Figures ||a and 26 show the X-ray radial surface 
brightness profiles and the best fit composite models for 
MS 0451 and CI 0016, respectively. The models for both 
clusters show a good fit to the data over a large range of 
angular radii. Using equation with the best-fit param- 
eters from Table || and the cooling functions from Table 0, 
we find the distance to MS 0451 to be 1278± 2 ^ Mpc and 



the distance to CI 0016 to be 20411^ Mpc, where the 
uncertainties are statistical only (see discussion below). 

Our fitting results are consistent with previous analyses 
of the ROSAT data of MS 0451 and CI 0016. Donahue 

(1996) analyzed the MS 0451 PSPC data and found j3 = 
1.01±^i| and e c = 38.2i x ^ arcseconds. Table | shows 
the comparison for CI 0016 with Neumann & Bohringer 

(1997) and Hughes & Birkinshaw (1998). 

We also compare our results with the distance determi- 
nation to CI 0016 by Hughes & Birkinshaw (1998). They 
analyzed the PSPC observations to determine the ICM 
shape parameters and then used that model to extract the 
SZE central decrement from observations taken with the 
OVRO single dish 40 m telescope at 20.3 GHz. They ob- 
served seven points in a north-south scan through CI 0016. 
Beam switching was done using the 40 m dual-beam sys- 
tem, which provides two 1'.78 FWHM beams separated by 
7'. 15 in azimuth. The central decrement extracted from 
such scans depends on the adopted center of the SZE sig- 
nal as well as the adopted ICM shape parameters, (3 and 
9 C . Interferometric observations provide two-dimensional 
imaging information with accurate astrometry and there- 
fore provide information about the cluster center and the 
ICM shape parameters. We find a central decrement of 
— 1242+^Qg /iK remarkably consistent with the Hughes & 

Birkinshaw value of — 1201+ifg /xK converted to thermo- 
dynamic temperature at 30 GHz. Hughes & Birkinshaw 
found the distance to CI 0016 to be 1 863^549 Mpc in good 
agreement with ours, where the uncertainty is statistical 
only and we have corrected for the frequency dependence 
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Table 4 

Comparison of Cl 0016 Analyses 



9 C 

Reference Instrument (3 (arcsec) 



NB97 a PSPC 0.80 ™ 50.5l£n 

NB97 a HRI 0.68 ±°-% 38.5±|;g 

HB98 b PSPC Q-728±gig3i 40.7±|;1 

This work joint PSPC & HRI 0.749±^g 42.3l2'o 



a Neumann fc Bohringcr 1997 
b Hughes fc Birkinshaw 199£ 



of the SZE (f( x \ —2) and relativistic corrections. Our 
PSPC central surface brightness and cooling function are 
both lower than the Hughes & Birkinshaw values. This 
difference arises entirely from using a different bandpass 
for the analysis (0.5 — 2.0 keV versus 0.4 — 2.4 keV). How- 
ever, only the ratio of the surface brightness and cooling 
function enters into the distance calculation. Our ratio 
times [j,h/ Me (A eH versus A ee ) is 1.65 xlO 11 arcmin~ 2 cm~ 5 
in fortuitously good agreement with theirs, 1.64 x 10 11 
arcmin cm" 5 . 

There is a known correlation between the (3 and 8 C pa- 
rameters of the (3 model. One might think this correlation 
would make determinations of D A imprecise because D A is 
calculated from these very shape parameters of the ICM. 
Figure || illustrates this correlation and its effect on D A 
for MS 0451. The filled contours are the 1, 2, and 3 a 
AS confidence regions for [3 and 8 C jointly. The lines are 
contours of constant D A in Mpc. With our intcrferometric 
SZE data, the contours of constant D A lie roughly parallel 




35 40 

6 (arcsec) 



Fig. 3. — Confidence regions from the joint SZE and X-ray fit 
for MS 0451. The filled regions are 1, 2, and 3 sigma confidence 
regions for /3 and C jointly (AS = 2.3, 6.2, 11.8) and the cross 
marks the best fit /3 and 8 C . Solid lines are contours of angular 
diameter distance in Mpc. The Da contours lie roughly parallel to 
the f3-0c correlation, minimizing the effect of this correlation on the 
uncertainties of D A . 



to the (3-Q c correlation, minimizing the effect of this corre- 
lation on the uncertainties of D A . Figure ^ shows similar 
behavior for Cl 0016. Different observing techniques will 
result in different behavior. Contours of constant D A have 
been found to be roughly orthogonal to th e (3-9 c correla- 
tion for some single dish SZE observa tions ( Birkinshaw & 
Hughes 1994] |Birkinshaw et al. 199l|). 



Uncertainties in the angular diameter distance from the 
fit parameters are calculated by gridding in the interest- 
ing parameters to explore the AS likelihood space. The 
most important parameters in this calculation are AT , 
S x0 , (3, and C . Radio point sources and the cosmic X-ray 
background affect AT and S x0 , respectively. Therefore 
we grid in AT , S xo , f3, and 9 C allowing the X-ray back- 
grounds for the PSPC and HRI to float independently. To 
estimate the effect of the radio point sources, we find the 
best-fit parameter values with the point source flux fixed 
at its best fit value. We also run the grids for the point 
source flux fixed at the ±1 a values. The point sources in 
the MS 0451 and Cl 0016 fields are both far enough from 
the cluster center so that their flux contributions have a 
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Table 5 

D A Observational Uncertainty Budget (percent) 



Cluster 


Fit a 


N H h 


[Fe]/[H] C 


rp b 

e 


Total d 


MS 0451 
CI 0016 


+ 13.7 
-13.1 
+ 17.8 
-16.4 


+0.9 
-1.2 
+ 1.1 
-1.2 


+1.1 
-i.i 

+2.1 
-1.3 


+ 15.4 
-19.2 
+ 15.4 
-19.1 


+20.7 
-23.3 
+23.7 
-25.2 



a The 68.3% uncertainties over the four-dimensional 
error surface for f3, 9 C , S x0 , and AT . 

h D A decreases as parameter increases. 

c Metallicity relative to solar. 

d Combined in quadrature. 




45 50 

(arcsec) 

n 

Fig. 4.— Same as Figure |] for CI 0016. 

negligible effect on the central decrement and do not 
change the cluster shape parameters significantly. From 
this four dimensional AS hyper-surface, we construct con- 
fidence intervals for each parameter individually as well 
as confidence intervals for D A due to S xo , AT , (3, and 
6 C jointly. The correlations between the /3 model param- 
eters require this treatment to determine accurately the 
uncertainty in D A from the fitted parameters. To com- 
pute the 68.3% confidence region we find the minimum 
and maximum values of the parameter within a AS* of 1.0. 
We emphasize that these uncertainties are meaningful only 
within the context of the spherical isothermal (5 model. 

The observational uncertainty budget for D A is shown 
in Table [s]. The uncertainties in the fitted parameters 
come from the above procedure. The only other parame- 
ter that enters directly into the D A calculation is T eo . Since 
D A oc T eo ~ 2 , the uncertainty in D A due to T eo is listed as 
twice the fractional uncertainty on T e0 . The other param- 
eters, column density and metallicity, as well as T e0 , affect 
the X-ray cooling function. We estimate the uncertain- 



ties in D A due to these parameters by taking their 68.3% 
ranges and seeing how much they affect the cooling func- 
tion. The uncertainty in the cooling function due to T eo 
is <J 0.5% and is ignored. The uncertainty on D A due to 
observations is dominated by the uncertainty in the elec- 
tron temperature and the SZE central decrement. Note 
that changes of factors of two in metallicity result in a 
~ 1% effect on D A . The column densities measured from 
the X-ray spectra are differ ent from those from H I surveys 
(Dickey & Lockman 1990). We use the column densities 
from X-ray spectral fits since that includes contributions 
from nonneutral hydrogen and other elements which ab- 
sorb X-rays. The survey derived column densities change 
the angular diameter distance by ~ ±5%, which we in- 
clude as a systematic uncertainty (see § ||). 

To determine the Hubble Constant, we perform a y 2 
fit to our calculated D A s versus z for three different cos- 
mologies. To estimate statistical uncertainties, we com- 
bine the uncertainties on D A listed in Table ^ in quadra- 
ture, which is only strictly valid for Gaussian distributions. 
This combined statistical uncertainty is symmetrized (av- 
eraged) and used in the fit. We find 



H = I 57 



52t J ° km s- 1 Mpc" 1 
km s _1 Mpc -1 
Mpc -1 



+11 

- 8 
j + 12 



km s- 1 



n M =i-o, O\=o.o, 

fi M =0.3, tt A =0.0, (7) 
n M =0.3, O\=0.7, 



where the uncertainties are statistical only. The statistical 
error comes from the \ 2 analysis and includes uncertain- 
ties from T e , the parameter fitting, metallicity, and N H (see 
Table ||) . We have chosen three cosmologies encompassing 
the currently favored models. There is a ~ 20% range in 
H at z ~ 0.5 due to the geometry of the universe. 

4.1. Sources of Possible Systematic Uncertainty 

The absolute calibration of both the SZE observations 
and the PSPC and HRI directly affects the distance deter- 
minations. The absolute calibration of the interferomet- 
ric observations is conservatively known to about 4% at 
68% confidence, corresponding to a 8% uncertainty in H 
(oc A77 2 ). The effective areas of the PSPC and HRI are 
thought to be known to about 10%, introducing a 10% 
uncertainty into the H Q determination through the calcu- 



10 



Table 6 

H Systematic Uncertainty Budget for MS 0451 and Cl 0016 



Systematic 


Eitcct 




(%) 


SZE calibration 


± 8 


X-ray calibration 


±10 


N H " 


± 5 


Asphericity a 


±14 


Isothermality & clumping 


±20 


Undetected radio sources' 3 


±16 


Kinetic SZE a 


± 6 


Total 


±33 



a Includes a l/v2 factor for our 2 
cluster sample. 

b Average of effect from the two 
cluster fields. 

c Combincd in quadrature. 



lation of S. In addition to the absolute calibration uncer- 
tainty from the observations, there are possible sources of 
systematic uncertainty that depend on the physical state 
of the ICM and other sources that can contaminate the 
cluster SZE emission. Table ^ summarizes the systematic 
uncertainties in the Hubble constant determined from MS 
0451 and Cl 0016. 

4.1.1. Cluster Atmospheres and Morphology 

Most clusters do not appear circular in radio, X-rays, or 
optical. Fitting a projected elliptical isothermal (3 model 
gives an axial ratio of — 0.80 and — 0.84 for MS 0451 and 
Cl 0016, respectiv ely, close to the local average of 0.80 
(Mohr et al. 1995). Under the assumption of axisymetric 
clusters, the combined effect of cluster asphericity and its 
orientation on the sky cons pires to introduce a ~ ±20% 
random uncertainty in H (Hughes fe Birkinshaw 1998|). 



When one considers a large, unbiased sample of clusters, 
presumably with random orientations, the uncertainty due 
to imposing a spherical model will cancel, manifesting it- 
self in the statistical uncertainty and allowing a precise de- 
termination of H . Recently, Sulkanen (1999) studied pro- 
jection effects using triaxial (3 models. Fitting these with 
spherical models he found that the Hubble constant esti- 
mated from the fitting was within ~ 5% of the true value. 
We are in the process of using N-body and smoothed par- 
ticle hydrodynamics (SPH) simulations of 48 clusters to 
quantify the effects of complex cluster structure on our 
results. 

Cooling flows also affect the derived distance to the clus- 
ter, affecting the emission weighted mean temperature and 
enhancing the X-r ay central surface brightness (see, e.g., 
Nagai et al. 2000|). A characteristic cooling time for the 



KIM 



is 



cool 



3kT e n tot 3kT e fi H 



(8) 

2An e n H 2An e fi tot ' 
where A is the bolometric cooling function of the cluster 



and all quantities are evaluated at the center of the clus- 
ter. Cooling flows may occur if the cooling time is less than 
the age of the cluster, which we conservatively estimate to 
be the age of the universe at the redshift of observation, 
t C ooi < t}i{z). For a flat, Einstein-de Sitter universe, the 
Hubble time is t H (z) = |iJ _1 (1 + z) _3/2 . Both MS 0451 
and Cl 0016 are observed in the z ss 0.55 universe so that 
tjj(z = 0.55) = 3 x 10 9 years. The ratio of the cooling 
time to the Hubble time for typical ICM parameters at 
redshift of 0.55 is then 



t H (z = 0.55) 



I9h 



; keV/ 



2.5 x 10 23 erg cm 3 s 1 



A 



io- 



(9) 



Using the best fit parameters we find A = 2.8 x 10 -23 
erg cm 3 s -1 and n e = 1.3 x 10~ 2 cm~ 3 for MS 0451 and 
A = 2.4 x 10~ 23 erg cm 3 s _1 and n e = 7.0 x 10~ 3 cm -3 for 
Cl 0016 (the densities are determined by eliminating D A 
in eqs. M and H in favor of n eo ). This implies t coo i/tu 
ratios of ~ 1.7 and ~ 2.6 respectively. These ratios are 
summarized in Table ^ for all three cosmologies considered 
in this paper. From this simple calculation, we do not ex- 
pect cooling flows in cither of these clusters. The X-ray 
radial surface brightness profiles (Figs ||a and b) provide 
no evidence for excess emission in the cluster core (see also 
Donahue fc Stocke 1995| ;|N cumann fc Bohringer 1997| ). As 
a check, we calculate t coo i/t[i ratios for each cluster ana- 
lyzed by Mohr et al. (1999). We check our cooling flow 
and non-cooling flow determinations versus those of Peres 
et al. (1998) and Fabian (1994). Of the 45 clusters in the 
Mohr sample, 41 have published mass deposition rates. 
We assume the cluster does not contain a cooling flow if 
its mass deposition rate is consistent with zero, otherwise 
it is designated as a cooling flow cluster. We are able to 
predict whether a cluster has a cooling flow or not with a 
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90% success rate, suggesting that the ratio t coo i/tH pre- 
sented in equation (m is a good predictor for the presence 
of a cooling flow. 

An isothermal analysis of a non-isothermal cluster could 
result in a large distance error; moreover, an isothermal 
analysis of a large cluster sample could lead to systematic 
errors in the derived Hubble parameter if most clusters 
have similar departures from isothcrmality ( Birkinshaw & 
Hugh es 1994| [Inagaki et al. 1995j |Holzapfcl ct al. 1997b| ). 
The effects of temperature variations depend on the ob- 
serving technique. For example, PSPC X-ray constraints 
on the ICM distribution are very insensitive to tempera- 
ture variations for gas whose tempera tures are above 1.5 
keV (see Fig. 1, |Mathiesen ct al. 1999[ ). In principle, SZE 
observations are sensitive to temperature variations, be- 
cause the SZE decrement is proportional to the projected 
pressure distribution (see eq. 0). However, interferomet- 
ric observations of the type presented here are relatively 
insensitive to modest ICM temperature variations. 

Interferometric SZE observations sample the Fourier 
transform of the sky brightness distribution over a limited 
region of the u-v plane. Specifically, the radio telescope 
dish size imposes a minimum separation for any two tele- 
scopes, making it impossible to sample the Fourier trans- 
form of the cluster SZE below some minimum radius in the 
u-v plane. Moreover, the primary beam of the telescope 
defines some effective field of view, making interferomet- 
ric observations completely insensitive to sky brightness 
fluctuations on any angular scale for those regions of the 
sky which lie outside the field of view. For these reasons 
interferometric SZE observations are insensitive to large 
angular scale variations in sky brightness. Therefore, clus- 
ters whose core regions are approximately isothermal and 
whose ICM temperatures decrease only gradually toward 
the virial region pose no problems for an isothermal anal- 
ysis. 

We are currently analyzing mock observations of gas- 
dynamical cluster simulations to explore the effects of ex- 
pected temperature distributions for our observing strat- 
egy. These simulated clusters exhibit X-ray merger signa- 
tures consistent with those observed in real clusters and, 
presumably, they exhibit the appropriate complexities in 
their temperature structure as well. Preliminary results 
from this analysis indicate that expected temperature gra- 
dients do not introduce a large systematic error in our dis- 
tance measurements. However, clumping within the ICM 
due to the common mergers of subclusters does enhance 
the X-ray surface brightness by ~ 20%. This enhance- 
ment cau ses X-ray gas mass estimates to be biased high 
by 10% ( [Mohr et al. 1999b| ), and it results in a ~ 20% 
underestimate of cluster distances. There is currently no 
direct observational evidence of clu mping within the ICM, 
but merger signatures are common (Mohr et al. 1995), and 
the mergers are the driving mecha nism behind these fluc - 
tuations in the simulated clusters ( Mathicscn ct al. 1999| ). 
We conservatively include a 20% systematic for clumping 
and departures from isothermality. 

4.1.2. Possible SZE Contaminants 

Undetected point sources may bias the angular diameter 
distance. Point sources near the cluster center mask the 
SZE decrement, causing an underestimate in the magni- 
tude of the decrement, and therefore an underestimate of 



the angular diameter distance. While we can not rule out 
point sources below our detection threshold, to estimate 
an upper bound on their effects we add a point source 
with flux at our detection limit near the cluster center 
to each data set and then fit the new data set, not ac- 
counting for the added point source. Such a point source 
being at the cluster center is highly unlikely but provides 
an upper bound to the effects of undetected point sources. 
For a point source with flux density 1, 2, and 3 times 
the rms (~ 65 ^Jy) in the high-resolution (> 1,800A) 
image for MS 0451, we find the magnitude of the decre- 
ment decreases by 3%, 8%, and 14% respectively. For 
the CI 0016 high-resolution image (;> 1, 500A), we find the 
decrement changes by 2%, 10%, and 17% for an on-center 
point source with flux density 1, 2, and 3 times the rms 
(~ 90 ^tJy) respectively. However, we can place more strin- 
gent constraints on contamination from undetected point 
sources because we have information about the distribu- 
tion of point sources in these two fields from observations 
at lower frequencies. By performing a deeper survey of 
these fields with our 30 GHz receivers with an array con- 
figured for higher resolution we can lower our point source 
detection threshold until the uncertainty from undetected 
point sources b ecomes negligible. 

The NVSS (|Condon ct al. 1998|) detected two point 
sources within 400" from the center of MS 0451. Of these, 
we detect the one that is 170" from the pointing center, 
but the point source 295" from the pointing center is out- 
side the OVRO field of view. Sources with flux densities 
greater than ~ 4 a ss 2 mJy appear in their catalog. As a 
more realistic upper bound on the contamination from un- 
detected point sources, we extrapolate a point source with 
flux density equal to 4 a at 1.4 GHz to 28.5 GHz using 
the average sp ectral index of radio sources in galaxy clus- 
ters a = 0.77 ( Cooray et al. 1998 ). Therefore, we place a 
180 //Jy point source near the center of MS 0451 and then 
fit the new image, not accounting for the additional point 
source. The magnitude of the central decrement decreases 
by 13% which is a reasonable upper bound to the contami- 
nation from undetected point sources in the MS 0451 field 
and similar to the constraints derived from our own data. 

Moffet & Birkinshaw (1989) surveyed the region around 
CI 0016 at 5 GHz with the VLA and then followed up 
the 5 GHz sources at 1.4 GHz and 15 GHz. Three of their 
sources (10, 14, and 15) fall within the BIMA field of view. 
We detect source 15 in the BIMA data, but it falls outside 
the OVRO field of view at 338" from the pointing center. 
We extrapolate sources 10 and 14 to 28.5 GHz from the 
1.4 GHz observations using the spectral index a — 0.77, 
which is consistent with the Moffet & Birkinshaw result, 
a = 0.7. After correction for the primary beam, sources 
10 and 14 are expected to be 28 fiJy and 227 [iJy, respec- 
tively. We add these two point sources to the OVRO data 
placing them at their NVSS positions, perform a model fit 
not accounting for them, and find a 3% change in the cen- 
tral decrement. Moffet & Birkinshaw searched for peaks 
that were 5er or greater. Extrapolating the 5 GHz rms of 
80 /iJy to 28.5 GHz results in a 21 /j,Jy rms. Placing a 5<r 
(100 /xJy) point source near the cluster center decreases 
the magnitude of the central decrement by 3%, identical 
to the combined effect from sources 10 and 14. 

Cluster peculiar velocities with respect to the CMB in- 
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Table 7 
Ratio of t cool /t H (z) 





Cosmology (Om, 


n A ) 


Cluster 


(1.0, 0.0) (0.3, 0.0) 


(0.3, 0.7) 


MS 0451 


1.7 1.3 


1.0 


CI 0016 


2.6 1.9 


1.6 



troduce an additional CMB spectral distortion known as 
the kinetic SZE. The kinetic SZE is proportional to the 
thermal effect but has a different spectral signature so it 
can be disentangled from the thermal SZE with spectral 
SZE observations. For a 10 keV cluster with a line-of-sight 
peculiar velocity of 1000 km s _1 , the kinetic SZE is ~ 11% 
of the thermal SZE at 30 GHz. Watkins (1997) presented 
observational evidence suggesting a one-dimensional rms 
peculiar velocity of ~ 300 km s _1 f or clusters, and re 
cent simulations found similar results ( Colberg et al. 2000| ) 



With a line-of-sight peculiar velocity of 300 km s 
more typical 8 keV cluster, the kinetic SZE is 



1 and a 
4% of the 



thermal effect, introducing up to a ~ ±8% correction to 



the angular diameter distance computed from one cluster. 
The effects from peculiar velocities when averaged over an 
ensemble of clusters should cancel, manifesting itself as an 
additional statistical uncertainty similar to the effects of 
asphericity. 

CMB primary anisotropies have the same spectral sig- 
nature as the kinetic SZE. Recent BIMA observations 
provide limits on primary anisot ropics on the scales o f 
the observations presented here ( Holzapfel et al. 2000p . 
We place a 95% confidence upper limit to the primary 
CMB anisotropies of AT < 22 at I ~ 5500 (~ 2' 
scales). Thus primary CMB anisotropies are an unimpor- 
tant (<; 2%) source of uncertainty for our observations. 



and a 6% (~ 8/v2) uncertainty from the kinetic SZE. 
These systematic uncertainties are summarized in Table [| 
The uncertainty from undetected radio sources is the aver- 
age of the maximum effects due to undetected sources for 
the MS 0451 (26%) and CI 0016 (6%) fields. The contri- 
butions from asphericity and kinetic SZE should average 
out for a large sample. 

Our H a determination from high-redshift clusters is 
consistent with other SZE based measurements as well 
as the recent results from the Hubble Space Telescope 
(HST) H Key Project, which probed the near by uni- 
vcrse and f ound H Q = 71 ± 6 km s^ 1 Mpc" 1 QMould 
et al. 2000"|). Birkinshaw (1999) compiled current SZE 



based H measurements and found an ensemble average 
of ~ 60 km s _1 Mpc -1 independent of the chosen cosmol- 
ogy; where the uncertainty is difficult to ascertain because 
the measurements are not independent, many share SZE 
or X-ray data, and nearly all share common absolute cal- 
ibrations. 

The SZE derived distances are direct, making them an 
interesting check of the cosmological distance ladder. Re- 
cent observations of mas ers orbiting the nucle us of the 



nearby galaxy NGC4258 (Herrnstein et al. 1999) illustrate 



a method of determining direct distances in the nearby 
universe. Time delays from analysis of gravitational lens- 
ing data from galaxy clusters are another direct distance 



5. DISCUSSION AND CONCLUSIONS 

We perform a maximum-likelihood joint fit to interfer- 
ometric SZE and ROSAT X-ray (PSPC and HRI) data 
to constrain the ICM parameters for MS 0451 and CI 
0016. We model the ICM as a spherical, isothermal f3 
modcL — From this analysis we determine the distances 



to KpI197«+265 Mpr. and 2041+^ Mpr fnr MS fUM artrl 

CI 0016, respectively (statistical uncertainties only). To- 
gether, these distances imply a Hubble constant of 



52 



km s" 1 Mpc- 1 ; Q M =l-0, ft A =0.0, 
H = { 57t 1 g t\l km s- 1 Mpc" 1 ; tl M =0.3, f2 A =0.0, 



-10 +17 
- 7 -17 



63 



8 -19 
-12 +21 
- 9 -2 



[ km s- 1 Mpc- 1 ; Q M =0.3, ft A =0.7, 

( 10 ) 

where the uncertainties are statistical followed by system- 
atic at 68% confidence. The systematic uncertainties have 
been added in quadrature and include an 8% (4% in AT ) 
uncertainty from the absolute calibration of the SZE data, 
a 10% effective area uncertainty for the PSPC and HRI, a 
5% uncertainty from the column density, a 14% (~ 20/ v2) 
uncertainty due to asphericity, a 20% effect for our as- 
sumptions of isothermality and single-phase g 16% 
% in AT ) uncertainty from undetected radio sources, 



cent examples see 


Fassnacht ct al. 1999; 


Biggs ct al. 1999 


Lovell et al. 1998|; parkana 1997; |Schec 


titer et al. 1997) 



The redshift independence of the SZE makes it a power 
ful probe of clusters at high redshift. The combination of 
SZE and deep X-ray observations could provide a valuable 
independen t check of high redshift S N la results ( [Schmidt 



et al. 1998; Pcrlmutter et al. 1999), which constrain the 



geometry of the universe. 

We are currently analyzing a larger sample of SZE clus- 
ters which will reduce the statistical uncertainty as well 
as effects from asphericity and the kinetic SZE. Analysis 
of mock observations of simulated clusters will also pro- 
vide insight into the effects of temperature gradients and 
multiphase ICMs. With the recent launch of Chandra and 
the impending launch of XMM, we will soon obtain bet- 
ter T e measurements (currently a large source of observa- 
tional uncertainty; see Table ra) and measure temperature 
profiles. In addition, the ~ 2% absolute calibration un- 
certainty of Chandra will soon replace the ~ 10% ROSAT 
absolute calibration uncertainty. There is also work being 
done to improve the absolute calibration at 30 GHz using 
the planets. The goal is to achieve a <; 1% absolute cal- 
ibration, further reducing the systematic uncertainties in 
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the derived Hubble parameter. 
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